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ABSTRACT 

The first stars are the key to the formation of primitive galaxies, early cosmological reionization 
and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their 
extreme luminosities, individual Population III stars will likely remain beyond the reach of direct 
observation for decades to come. However, their properties could be revealed by their supernova 
explosions, which may soon be detected by a new generation of NIR observatories such as JWST 
and WFIRST. We present light curves and spectra for Pop HI pair-instability supernovae calculated 
with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for 
the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and 
Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR 
signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae 
out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan- 
STARRS will find them at z < 7 — 8. The discovery of these ancient explosions will probe the first 
stellar populations and reveal the existence of primitive galaxies that might not otherwise have been 
detected. 

Subject headings: early universe - galaxies: high-redshift - stars: early-type - supernovae: general - 
radiative transfer - hydrodynamics - shocks 



1. INTRODUCTION 

Population HI (Pop HI) stars are the key to under- 
standing primeval galaxies, the chemical enrichment and 
reionization of the early IGM, and the origin of super- 
massive black holes. Unfortunately, because these stars 
are at the edge of the universe there are no observational 
constraints on their masse s. Even with their extreme lu- 
minosities (jSchaereii |2002|) and the advent of new near 
infrared (NIR) observatories such as the James Webb 
Space Telescope {JWST), the Wide-Field Infrared Sur- 
vey Telescope ( WFIRST) and the Thirty-Meter Tele- 
scope (TMT), individual Pop III stars will remain be- 
yond the reach of o bservation for decades to come (e.g. 
[Gardner et"al]|2006l but see Rydberg ct al. 2012 for the 
possibility of detection by strong gravitational lensing). 

Numerical simulations indicate that primordial stars 
are born in 10^ - 10^ Mq dark matter halos at 
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z ~ 20 - 30. The original calculations suggested 
that they are 30 - 500 M(7^ and for m in iso l ation , 
one per halo (iBromm et all [TQOQ: Abel et all 1 2000 . 



20021: iBromm e t al. 2002: Naka mura fc Umemural 12001 . 
O'Shea fc Norman ,2007. : .Yoshida et all l2008f ). Newer 



models hav e since shown tha t some Pop HI stars form 
in binaries (jTurk et al.l 120091 ) . and that some may form 
with much lower masses in small clusters of up to a 
doze n (IStacv et all [20T0I : IClark et al.l IMlt IGreif et al.l 
I2011ll2012f) . Nevertheless, in spite of their increasing so- 
phistication these simulations fall well short of predict- 
ing the Pop HI initial mass function (IMF) because none 
of them realistically bridge the gap between the forma- 
tion and fragmentation of the protoste llar disk and its 
photoevapq r ation up to a Myr later (lOmukai fc Pallal 
2OOI 12001 iTan fc McKed [2003: IMcKee fc TanI l200a 



Stacy et alll2012l: iHosokawa et aLir2011| ). 

Some have tried to estimate the masses of Pop HI stars 
by reconciling their nucleosynthetic yields with the chem- 
ical abundances found in ancient, dim metal-poor stars 
in the Galactic halo, some of which may have formed 
in the ashes of the first stars (Beers fc Christliebl 120051 : 
iFrebel et all [20051 : iCaffau et aLll2012[ ). The current con- 
sensus is that 15-40 Mq Pop III stars die in core-collapse 
supernovae (CC SNe) and 140 - 260 Mq stars explode as 
far more energetic pair- instability supernovae (PI SNe), 
with 1 00 times the energy of Type la or Type II explo- 
sions (jHeger fc Wooslevll2003) . The metals expelled by 
the first SNe into the IGM depend on how mixing and 
fallback onto compact remnants occurs inside the star 
rior to shock breakout from its surface. iJoggerst et al.l 
2OT0t ) (hereafter JETIO) find that in 15 - 40 Mq Pop III 
SNe the shock fully disrupts the shells of elements in the 
star, heavily mixing the ejecta by the time it ruptures 
the surface. In PI SNe the shells expand almost homol- 
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ogously, with only occasional minor mixing between the 
O an d He laye rs prior to breako ut (iJoggerst fc WhalenI 
[20H (see also lChen et"aI|[20Tl[) . 

JET 10 find that the aggregate yields of 15 - 40 Mq Pop 
III CC SNc arc consistent with the chemical abundances 
measured in a. sample of ^ 130 extremely meta l-poor 
(EMP) stars (|Cavrel et al.|[200l iLai et alj [20081 ) . sug- 
gesting that some primordial stars may have much lower 
masses than inferred in the first numerical simulations. 
These results together with the failure to detect the 
odd-even nucleosynthetic pattern of PI SNe in extremely 
metal poor stars to date have led some to conclude that 
Pop III were only tens of solar masses, not hundreds, 
and that Pop III PI SNe may not have occurred in the 
early universe. However, such conclusions are premature 
for several reasons. First, the gas in which later genera- 
tions of stars formed may have been enriched by PI SNc 
to metallicities above those tar geted by surveys thus far 
(jKarlsson et all 120081) (see also lRen et all 120121 . for de- 
tails on a new PI SN survey). Second, extracting the 
Pop HI IMF from the fossil abundance record remains 
problematic. Larger samples of EMP stars are needed 
to better constrain SN yields, systematic error has been 
found in measurements of some elements, and how the 
first metals are taken up into new stars is not fully un- 
derstood. The recent discovery that rotating Pop HI 
stars down to 65 Mq may die as PI SNe significantly 
increases the likelihood of such events at high redshift 
(jChatzopoulos fc Wheeled l2012f ). Finally, and perhaps 
most compellingly, P I SNe have been discovered today 
(|Gal-Yam et al.ll2009[ ) in environments that are far less 
favorable to the for mation of massive progenitors than in 
the early universe (IPan et al.ll20T2B . 

The properties of the first stars may soon be unveiled 
by their SN explosions, which could be found by JWST 
or WFIRST because they are thousands of times brighter 
than their progenitors or their host galaxies. Past stud- 
ies have considered detection thresholds for PI SNe fo r 
various instrument s at z ^ 5 (fScannapieco et al.ll2005D . 
6 < z < 15 (P an et al.ll2012aD . and in rough terms out 
to z - 30 (iHummel et al.ll2012| ). However, to accurately 
model the NIR signatures of these events, PI SNe must 
be simulated with radiation hydrodynamics in realistic 
circumstellar envelopes and their spectra must be con- 
volved with Lyman absorption by the neutral IGM at 
early epochs and redshifted. The interaction of the SN 
with the envelope shrouding the star and the opacity of 
the envelope itself both crucially shape the spectra in 
the source frame and govern how much fiux is eventually 
redshifted into the NIR. Lyman absorption by neutral 
hydrogen prior to the era of reionization determines how 
much flux from these ancient SNc is removed or scat- 
tered from our line of sight. Previous studies, includ- 
ing earlier calculat ions of PI SN light curves and spectra 
(jKasen et al.ll20HT) . did not incor porate these ef f ects. 

Accounting for these processes, iWhalen et all (l2012bD 
recently found that JWST wih detect PI SNe out to z ~ 
30 and WFIRST will detect them out to z - 15 - 20. 
However, their grid of models did not cover the full range 
of stellar structures expected for 140 - 260 M© Pop HI 
stars. We have now calculated the NIR signals of Pop HI 
PI SNe for both blue compact and red hypcrgiant progen- 
itors from 7 < z < 30 with the Los Alamos RAGE and 
SPECTRUM codes. In § 2 we review the PI explosion 



mechanism, the physics that governs the presupernova 
structure of the star, and our Kepler Pop HI PI progen- 
itor and SN models. We describe how these explosions 
are evolved in RAGE and post processed in SPECTRUM 
to obtain source frame light curves and spectra in § 3 and 
we examine blast proflles and spectra in § 4. In § 5 we 
present NIR light curves for Pop HI PI SNe in the ob- 
server frame and determine their detection limits as a 
function of redshift. In §6 we calculate PI SN detection 
rates for 5 < z < 30, and in § 7 we conclude. 

2. PI SN MODELS 

PI SNe were firs t prop osed bv lRakavv fc Shavivl ()1967f ) 
and iBarkat et al.l (I1967D and have been stud ied numer- 
ous times since (see iHeger fc WooslevI l2002l and refer- 
ences therein). Pop HI stars above 100 Mq encounter 
the pair instability after central carbon burning, when 
thermal energy creates e^e~ pairs rather than maintain- 
ing pressure support against collapse. The cores of these 
stars subsequently contract, triggering explosive ther- 
monuclear burning of O and Si. Above 140 AIq, the 
energy that is released completely unbinds the star, and 
no black hole is formed. At 260 Mq the core of the star 
reaches temperatures that are high enough for alpha par- 
ticles to be photo-disintegrated into free nucleons, which 
consumes as much energy per unit mass as was released 
by all preceding burning. The star collapses instead of 
exploding. PI SNe are the most energetic thermonuclear 
events in the universe, with yields of up to 10^'^ erg for 
260 Mq stars. 

The early spectral signatures of Pop HI PI SNe heavily 
depend on the radius and structure of the star in addi- 
tion to the explosion energy, the interaction of the blast 
with the envelope, the opacity of the envelope, and ab- 
sorption by the IGM. Shock temperatures at breakout 
are lower for large stars than for compact stars of equal 
mass because the shock has expanded to a greater radius 
and done more work on its environment. Early spectra 
from the explosions of compact blue stars are therefore 
harder than those of red giants. The structure of the star 
determines the degree to which its interior is mixed be- 
fore the shock breaks out from its surface and hence the 
lines that appear in its spectra. The size and structure 
of the progenitor in turn are governed by its metallicity 
and by internal convection over its life. 

2.1. Semi-Convective Mixing 

Convection can determine if a very massive Pop HI 
star dies as a compact b lue giant or a red hyperg iant. 
As described in detail in iScannapieco et al.l ()2005l) . the 
convection zone of the star can approach, touch or even 
penetrate the lower hydrogen layers, mixing them with 
carbon dredged up from the core from He burning. When 
these two high temperature components mix, they burn 
violently, boosting energy production rates in the H shell 
by up to several orders of magnitude. This, together with 
the now greater opacity of the lower hydrogen layer, can 
puff up the star by more than an order of magnitude in 
radius. 

2.2. Metallicity 

Gas in high-redshift halos that is enriched to metallic- 
ities below lO"'' ''' Zq fragments on mass scales that are 
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Fig. 1. — Kepler PI SN explosion and stellar envelope profiles. Left panel: densities for all nine models. The surfaee of each star is 
marked by the abrupt drop in density by ten orders of magnitude at r ~ 10^^ em in the four z-series profiles and a few X 10^* cm in the 
five u-series stars. A wind with a total mass of 0.1 Mq, a velocity of 10* cm/s, and a free-streaming profile is extended from the 
surface of the star out to the uniform relic H II region surrounding the star, which is not visible on this plot. Right panel: velocity profiles 
for the u-series stars. Both the outgoing shock and ongoing collapse of the outer layers of the star are visible. 



essentially identical to those of pristine gas and forms 
very massive stars (e.g.[Bro mm et al. 2001; Mackc v et al.l 
120031 iSmith fc Sigurdssoi]|l2007t ). However, these small 
metallicities are more than enough to enhance CNO reac- 
tion and energy generation rates in the hydrogen burning 
layers of the star, inflating the stellar envelope as much as 
convection. Since there is a strong degeneracy between 
the effects of metals and convection on the structure of 
the star, the full range of light curves and spectra for 
Pop III PI SNe can be as easily spanned by metallicity 
as convective overshoot, as we demonstrate below. 

2.3. Explosive Mixing 

Mixing within the star during the explosion can de- 
termine which elements appear in emission lines over 
time. In the frame of the shock, the photosphere from 
which photons escape descends into the ejecta over time 
as each fluid element expands and the ejecta is diluted. 
If Rayleigh- Taylor (RT) mixing precedes shock breakout, 
it can dredge heavy elements up from a deeper mass co- 
ordinate to the photosphere and expose them to the IGM 
at much earlier times. Since mixing is minimal in Pop 
III PI SNe, the absence of Fc or Ni lines soon after shock 
breakout would be one of several markers of the event. 
The evidence or absence of mixing in early spectra is one 
manner in which a rudimentary IMF could begin to be 
built up by detections of Pop III explosions. 

2.4. Kepler 

To capture the structure of the progenitor, we evolve 
150, 175, 200, 225, and 250 Mq zero-metallicity stars 
(z-series) and 10^** Zq stars (u-series) from the zero 
age main sequence to the onset of colla pse in the ID 
Lagrangian stellar evolutio n code Kepler (jWeaver et al.l 
[1971 iWooslev et al|[200l . The SN begins when this 
collapse triggers explosive O and Si burning. Unlike the 
CC SN simulations of JETIO, in which the blast must 



be artiflcially launched with a piston and the explosion 
energy is a free parameter, the PI SN is an emergent fea- 
ture of our stellar evolution model, and its energy is set 
by how much O and Si burns. The blast was followed 
until the end of all nuclear burning at ~ 20 s, when the 
shock was still deep inside the star. We calculate energy 
generation with a 19-isotope network up to the point 
of oxygen depletion in the core and with a 128-isotope 
quasi-equilibrium network thereafter. The z-series 150 
Mq star collapses to a black hole without an explosion. 
The number of mass zones on the grid ranged from 1000 - 
1200 and was always sufficient to resolve all salient struc- 
tures of the star and SN. We show density and velocity 
profiles for our explosions in Figure [T] and summarize our 
grid of models in Table [T] As in all previous studies of 
PI SNe, we consider only non-rotating progenitors. 

We use metallicity rather than convective overshoot to 
parametrize the structures of our progenitors because for 
any given star the size of the central convection zone is 
uncertain but Z = and lO^"' Zq bracket the metallic- 
ities over which 140 - 260 M q Pop III stars are thought 
to form (iMackev et al.l 120031: ISmith fc Sigurdsso^ 120071 : 
ISmith et all 120091) ■ As we show in Figure [TJ these two 
metallicities yield compact blue giant and red hypcrgiant 
envelopes similar to those obtained by varying scmi- 
convective mixing in [Scannapieco et al.„ (200S)- AH u- 
series stars in our study die as red hypergiants and all 
z-series stars die as compact blue giants. The light curves 
and spectra for our SNe should therefore bracket those 
that will be found in observations. However, we note that 
most 140 - 260 Mq Pop HI stars are expected to exhibit 
convective mixing and die as red stars rather than blue 
stars. 

3. RAGE AND SPECTRUM SIMULATIONS 

We propagate the shock through the interior of the 
star, its surface, and then out into the surrounding 
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TABLE 1 

Kepler PI SN Models (masses are in Mq) 



model 


MHe 


R (10^^ cm) 


E (10^^ erg) 




ul50 


72 


16.2 


9.0 


0.07 


ul75 


84.4 


17.4 


21.3 


0.70 


u200 


96.7 


18.4 


33 


5.09 


u225 


103.5 


33.3 


46.7 


16.5 


u250 


124 


22.5 


69.2 


37.9 


zl75 


84.3 


0.62 


14.6 





z200 


96.9 


0.66 


27.8 


1.9 


z225 


110.1 


0.98 


42.5 


8.73 


z250 


123.5 


1.31 


63.2 


23.1 



medium with the rad iation hydrodynamics code RAGE 
(|Gittings et al.ll2008[) . RAGE (Radiation Adaptive Grid 
Eulerian) is a multidimensional adaptive mesh refine- 
ment (AMR) radiation hydrodynamics code developed at 
Los Alamos National Laboratory (LANL). RAGE cou- 
ples second-order conservative Godunov hydrodynam- 
ics to grey or multigroup flux-limited diffusion (FLD) 
to model strongly radiating flows. RAGE utilizes the 
extensive LANL O PLIB database of atomic opacitie^ 
(|Magee et al.l I1995D and can also evolve multimaterial 
flows with a variety of equations of state (EOS). We de- 
scribe the physics implemented in our RAGE models and 
why it is nee ded to capture the features of our light curves 
in detail in iFrev et all ()2012D (hereafter FET12): mul- 
tispecies advection, grey FLD radiation transport with 
2-temperature (2T) physics, energy deposition from ra- 
dioactive decay of ^^Ni, and no self-gravity. In particu- 
lar, 2T radiation transport, in which radiation and mat- 
ter temperatures are evolved separately, better models 
shock breakout and its aftermath when matter and ra- 
diation can be out of equilibrium. This is an important 
improvement over earlier IT models of PI SN explosions. 
We evolve mass fractions for 15 elements, the even num- 
bered elements predominantly synthesized in PI SNe. 

3.1. Initial Grid 

We map densities, velocities, specific internal energies 
(erg gm~^), and species mass fractions for the explosion 
and star from Kepler together with a circumstellar en- 
velope onto a 50,000 zone ID spherical AMR mesh in 
RAGE. Since radiation energy densities were not explic- 
itly evolved in the Kepler models, we initialize them in 
RAGE as 

-aT\ 
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where a = 7.564 x 10~^^ erg cm~^ K~^ is the radiation 
constant and T is the gas temperature. Also, because the 
gas energy in Kepler includes contributions by ionization 
states of atoms, we unambiguously construct the specific 
internal energy from T with 



CvT, 



(2) 



where Cy = 7.919 x 10^ erg K~^ is the specific heat 
of the gas. Because there is little mixing in the star 
as the shock propagates to the surface, the radial dis- 
tribution of elements in the star is essentially frozen in 
mass coordinate at death and expands homologously into 

http:/ /aphysics2/www. t4.lanl.gov/cgi-bin/opacity/tops. pi 



space thereafter. One-dimensional simulations are there- 
fore sufficient to capture all key attributes of our PI SN 
light curves and spectra. 

At the beginning of the simulation we allocate 5000 
zones from the center of the grid to the edge of the shock 
in the velocity profile on the root grid. We permit up to 
five levels of refinement in both the initial interpolation of 
the profiles onto this grid and throughout the simulation. 
This guarantees that all important density and velocity 
features in the initial blast profile are resolved by at least 
10 zones on the root grid and can be subsampled by up to 
32 times more zones as needed. It also ensures that the 
photosphere of the shock is always fully resolved, since 
failure to do so can lead to underestimates of luminosity 
during post processing. The shock radii in our Kepler 
profiles vary from 1.5 - 3.0 x 10^*^ cm and our initial 
root grid has a resolution of 1.0 x 10^ cm, with an outer 
boundary at 5.0 x 10^^ cm. 

We impose refiecting and outflow boundary conditions 
on the fluid and radiation flows at the inner and outer 
boundaries of the mesh, respectively. When the calcu- 
lation is begun, Courant times are initially small due to 
high temperatures, large velocities and small cell sizes. 
To reduce execution times and to accommodate the ex- 
pansion of the ejecta, we periodically regrid the profiles 
onto a larger mesh as the explosion grows, as described 
in detail in FET12. Each time we regrid the blast we 
allocate 5000 zones out to either the edge of the shock 
(pre-breakout) or to the edge of the radiation front (post- 
breakout - we take the point at which the radiation tem- 
perature falls below a few tenths of an eV to be the edge 
of the front). The inner boundary is always cm and the 
outer boundary of the final, largest mesh in our models 
is 1.0 X 10^* cm. Evolving the explosion out to 3 yr 
requires about 30 days on 32 processors on LANL plat- 
forms. 

3.2. Circumstellar Winds 

Just as there are no observational constraints on Pop 
III stars, there are none on their ambient environ- 
ments either. Numerical simulations of I- front break- 
out from primordial halos show that Pop III stars usu- 
ally drive all the baryons from their halos in super- 
sonic ionized fl ows and die in low-density H II regions 



(IWhalen et al.ll20 04: Kitava ma et al.| 



2004: 1 Alvarez et al.l 



l2006HAbel eranf2007i: IWise fc Ab ell 12008 al) . These relic 
H II regions have relatively flat density profiles n ~ 0.1 
- 1 cm^'^ out to 100 - 200 pc. Pop III stars are not 
thought to lose much mass over their lifetimes because 
there a re no metals in their atmospheres to drive strong 
winds fKudritzki' '2000'; 'Ba raffe et afl [200lt IVink et al.l 
120011: [Krticka & Kubat 200^. However, we allow for 
the possibility that they had some winds for several 
reasons. First, mass loss from very massive stars is 
usually just parametrized by metallicity, not calculated 
from first principle s with radiation hyd rodynamics and 
other physics (e.g.. iMevnet et aLlll994[) . Extrapolating 
zero wind from no metallicity could miss mass loss by 
other means, such as helium opacity in extreme lumi- 
nosities, hydromagnetic flows, or pulsational ejections 
(|Heger fc WooslevI [200l . Second, mixing in rotating 
Pop III stars, which we do not consider here, can dredge 
metals up from the interior of the star to its su rface and 
drive winds later in its life (e.g., lEkstrom et al- i2008i) . 
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although the mass loss is found to be minor. Finally, 
there are no observations to rule out winds from zero- 
mctallicity stars. We therefore extend a low-mass wind 
profile from the surface of the star out to the relic H II 
region: 



Pw{r) 



(3) 



where m is the mass loss rate due to the wind and vw is 
the wind speed. The mass loss rate is 



tuSL 



(4) 



where Aitot and tMSL are the total mass loss and main se- 
quence lifetime of the star, respectively. We take Mtot = 
0.1 M0, Vw = 1000 km s^^, H and He mass fractions to 
be 76% and 24%, and T = 0.01 eV in aU our models. 
The radius at which the density of the wind is smoothly 
joined to the H II region, whose density we take to be 
0.1 cm^^, varies with the radius of each star, but is typ- 
ically 10^^ - 10^® cm. Density and velocity profiles for 
the wind are shown in Figure [T] (where the sharp drop in 
density at the edge of the star where breakout occurs is 
also evident). If the wind density falls to the H II region 
density before the edge of the grid it is smoothly joined 
to that density at that radius. 

3.3. Ionization of the Wind 

To determine if the progenitor ionizes the wind we use 
the ZEUS-MP code to model the outgoing ionization 
front from t he surface of the star (jWhalen &: NormanI 
[20M [2008al lblV We center a blue 175 Mq z-series star 
and its wind on a ID spherical mesh with 200 zones and 
inner and outer boundaries at 7.0 x 10"'^^ cm and 1.0 x 
10^^ cm, the surface of the star and the outer regions 
of the wind. To enhance grid resolution in the densest 
regions of the wind we use logarithmically ratioed zones 
where 

= 1.043. (5) 

We propagate of the I-front from the star into the wind 
with multifrequency UV transport, with 40 bins uni- 
formly partitioned in energy from 0.255 - 13.6 eV and 80 
logarithmically spaced bins from 13.6 - 90 eV. We take 
the spectrum of the star to be blackbody, normalized 
to ionizing photon em ission rat e s, sur face temperatures, 
and luminosities from iSchaereJ ()2002D . 

We find that the blue star easily ionizes the wind on 
timescales of ~ lO'' yr. The gas temperature is 40,000 
K near the surface of the star, where rapid ionizations 
and recombinations in the large fluxes and densities there 
lead to greater heating than at 10^^ cm, where tempera- 
tures are ^ 25,000 K. Because the least massive blue star 
in our study easily ionizes the wind and has the highest 
surface wind density, we conclude that the winds around 
all z-series progenitors in our simulations are ionized. 

Red u-series progenitors have surface temperatures 
that are too low to emit ionizing UV but still ionize their 
wind envelopes because they are blue for most of their 
lives. They become too cool to sustain ionizing flux only 
in their final few hundred kyr. Recombination times in 



the wind therefore determine its ionization state when 
the star dies, without the need for a transport calcula- 
tion: 

^ (6) 



Here, ne is the electron number density and a{T) is 
the recombination rate coefficient for hydrogen, which 
we take to be 2.59 x 10"" T-° '^^ s-\ where T4 is the 
temperature in units of 10^ K. With the ionized gas tem- 
peratures we found for the compact 175 Mq star and the 
densities of the modest winds in our study, recombina- 
tion times vary from 5-50 kyr in the vicinity of the 
star, ensuring that the wind is neutral when the star 
explodes. For simplicity, we take the wind to be neu- 
tral around all the stars in our study, so our luminosities 
at shock breakout for z-series stars will be lower limits 
because neutral envelopes allow less flux to escape un- 
til they arc fully ionized by the breakout pulse. We note 
that the ionized wi nd itself has significan t luminosity due 
to recombinations (jRvdberg et al.l[2010( ) , but we arc pri- 
marily interested in the transient flux from the explosion, 
not emission from its envelope. 

3.4. SPECTRUM 

To calculate a spectrum from a RAGE profile, we map 
its densities, temperatures, velocities and mass fractions 
onto a new grid in the LANL SPECTRUM code. SPEC- 
TRUM performs a direct sum of the luminosity of every 
fluid element in this discretized snapshot to compute the 
total flux that exits the ejecta along the line of sight at 
every wavelength. This procedure is described in detail 
in FET12, and is formulated to account for Doppler shifts 
and time dilation due to the rclativistic expansion of the 
ejecta. We also calculate the intensities of emission lines 
and the attenuation of flux along the line of sight with 
detailed LANL OPLIB opacities, thereby capturing both 
limb darkening and the absorption lines imprinted on the 
flux by intervening material in the ejecta and wind. 

As explained in FET12, gas densities, velocities, mass 
fractions and radiation temperatures from the finest lev- 
els of refinement of the RAGE AMR grid are first ex- 
tracted and ordered by radius into separate files, one vari- 
able per file. These profiles typically contain more than 
50,000 radial zones, and constraints on machine memory 
and time prevent us from using all of them to calculate 
spectra, so only a subset of the points is mapped onto the 
SPECTRUM grid. We first sample the RAGE radiation 
energy density profile inward from the outer boundary 
to find the position of the radiation front, which we take 
to be where aT^ rises above 1.0 crg/cm"^. Wc then find 
the radius of the r = 20 surface by integrating the op- 
tical depth due to Thomson scattering inward from the 
outer boundary {hth = 0.288 for H and He gas at pri- 
mordial composition). This yields the greatest depth in 
the ejecta from which photons can escape because kth is 
the minimum opacity the photons would encounter. 

The extracted gas densities, velocities, temperatures 
and species mass fractions from RAGE are then interpo- 
lated onto a 2D grid in r and 9 in SPECTRUM whose 
inner and outer boundaries are and 10^^ cm, respec- 
tively. One hundred uniform zones in log radius arc allo- 
cated from the center of the grid to the r = 20 surface. 
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Fig. 2. — The highly radiative shock in the u225 model at 7.26 
X IC* s. The shock is at 1.2 X 10^** cm, deep inside the star. 
Photons from the shock cannot reach the surface of the star from 
this depth because of Thomson scattering in its upper layers. The 
0.3 eV plateau in the temperature at the stellar surface at 3.2 X 
10^* cm is optical radiation from the star propagating out into the 
wind before shock breakout. 



an intense burst of photons upon being exposed to the 
low density IGM, as shown in the right panel of Fig- 
ure [3l The advancing radiation front is visible as the 
plateau in gas temperature at 50 eV and 1.2 x 10"'^^ cm 
at 3990 s and 30 cV and 2.0 x lO^^ at 4550 s. The 
plateau temperatures are those to which the radiation 
front heats the ambient wind as it passes through it; the 
temperature of the shock that is emitting the radiation 
itself is much higher, as we show below. The temperature 
of the plateau falls as the shock expands, cools and its 
spectrum softens. Note that there are serious departures 
from self-similarity in the velocity profiles at late stages 
of breakout (third plot in the left panel of Figure ^ be- 
cause of significant coupling between radiation and gas 
in the shock. This is one reason Sedov- Ta ylor profiles 
are n ot good solutions for shock breakout (jFrver et al.l 
l2010f l and why radiation transport is required to model 
the flow. 

We show bolometric luminosities for the breakout tran- 
sient for the u-series and z-series PI SNe in the left and 
right panels of Figure|4l respectively. For comparison, we 
show the luminosity of the shock for the u225 explosion 
if it is assumed to be a blackbody. 



(7) 



The region between the r = 20 surface and the edge of 
the radiation front is partitioned into 5000 uniform zones 
in radius. The wind between the front and the outer 
boundary of the grid is then divided into 100 uniform 
zones in log radius, for a total of 5200 radial bins. The 
data in each of these new radial bins is mass-averaged 
to ensure that the SPECTRUM grid captures very sharp 
features from the original RAGE profile. The grid is uni- 
formly discretized into 160 bins in fi = cos0 from -1 to 
1. Our choice of mesh yielded good convergence in spec- 
trum tests, fully resolving regions of the flow from which 
photons can escape the ejecta and only lightly sampling 
those from which they cannot. 

Summing the luminosities at each wavelength in one 
spectrum yields the bolometric luminosity of the SN at 
that moment. Many such luminosities computed over 
a range of times constitutes the light curve of the ex- 
plosion. We sample shock breakout with 50 spectra uni- 
formly spaced in time and the rest of the light curve with 
another 200 spectra logarithmically distributed in time 
out to 3 yr. One SPECTRUM calculation averages 3-6 
hours on 32 processors on LANL machines. 

4. PI SN BLAST PROFILES, LIGHT CURVES AND 
SPECTRA 

The onset of explosive burning drives a strongly radi- 
ating shock from the O layer into the upper layers of the 
star, as we show in Figure [21 At this point the shock is 
not visible to an external observer because its photons 
are trapped by e~ scattering in the intervening layers, so 
they are simply advected outward by the fluid flow. 

4.1. Shock Breakout 

When the shock reaches the surface of the star it 
abruptly accelerates to large velocities in the steep den- 
sity gradient there, as we show in the left panel of Fig- 
ure [3] The acceleration heats the shock, and it releases 



where cr is the Stefan-Boltzmann constant and the tem- 
perature T of the shock is taken at the tth = 1 sur- 
face, where kth = 0.288. If photons of all energies broke 
through the surface of the star at the same time, the 
duration of the initial transient would be comparable to 
the light crossing time of the star, since photons emitted 
from its poles and its equator would reach an observer 
at times that differ by the time it takes light to cross the 
star. Given that u-series stars have radii of a few 10^^ 
cm, their breakout transients would last 1-2 hr, which 
is consistent with the width of the Airr^crT'^ approxima- 
tion to the luminosity of the pulse. In reality, the tran- 
sient is smeared out over longer times because photons 
at different energies escape the shock from a range of 
optical depths over time and because radiation remains 
strongly coupled to the shock past breakout (as shown in 
the third velocity profile in Figure[3]). The resulting pulse 
is dimmer but longer. The width of the transient is much 
greater for giant red u-scries stars because of their larger 
radii, lower surface densities, and thus broader regions 
from which photons break free of the shock. 

We show spectra for the breakout pulse for u250 and 
z250 PI SNe in the left and right panels of Figure [5l 
respectively. The z250 transient is mostly x-rays and 
the u250 pulse is both x-rays and hard UV. The spec- 
trum of the u250 transient is softer because the progen- 
itor is much larger, and the shock has done more PdV 
work on its surroundings and cooled more by the time it 
breaks out of the star. The ionization of the surrounding 
wind envelope is particularly evident in the z250 break- 
out spectra, in which prominent absorption features from 
100 to 10^ A in the first profile are mostly absent in the 
second. By the end of shock breakout the surrounding 
wind is completely ionized by the radiation pulse and the 
spectrum of the shock has essentially a blackbody pro- 
file. The absorption features are more prominent in the 
z250 spectra because of the higher wind densities at the 
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Fig. 3. — Shock breakout in the z225 explosion. Left panel; velocities, from left to right, at 3560 s, 4140 s, and 4550 s. Shock breakout 
through the surface of the star is evident in the abrupt jump in velocity in the second profile, and the continued acceleration of the shock 
down the density profile of the wind is visible in the third profile. Right panel; gas temperatures at 3560 s, 3990 s, and 4550 s from 
left to right. The radiation breakout pulse is the flat 50 eV plateau in gas temperature at 1.2 X 10^'^ cm at 4140 s. As the shock expands 
and cools the spectrum of the breakout pulse softens, which is why the gas temperature in its wake has fallen to ~ 30 eV at 4550 s. 
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surface of this smaller star. 

Peak bolometric luminosities vary from 7 x lO"'^ to 4 x 
10'''^ erg/s, or ~ 2500 times the luminosity of our galaxy. 
There is a general trend of greater luminosity with ex- 
plosion energy and later breakout times with progenitor 
mass, since stellar radii increase with mass. Although 
the shock has a much smaller radius at breakout in the z- 
scries than in the u-series, it has approximately the same 
total luminosity because it is hotter. This can be seen 
from how the luminosity i of a blackbody scales with 
radius and temperature in equation [71 Typical shock 
temperatures at the r = 1 surface at breakout are ^ 
1200 eV in the z-series stars and ^ 150 eV in the u-series 
stars, which roughly compensates for the factor of 10 in 



radius between the two progenitors. However, we note 
that the shock is not a blackbody and that its luminosity 
is roughly an order of magnitude below that predicted by 
equation[7l as we show in Figure IIFI . Treating SN shocks 
as blackbodics in general can lead to overestimates in lu- 
minosity of an order of magnitud e or more at breakout, 
as discussed in lFrver et al.l ( 20101 ). 

Wc also note that our transport solutions yield peak 
bolometric luminosit ies that are an ord er of magnitude 
higher than those of iKasen et al.l ()2011[ ). There may be 

Because of e~ scattering in the photosphere of the ejecta, the 
shock is more like a grey body whose luminosity L = Aner^aT'^ , 
where e < 1 is the correction to the blackbody luminosity due to 
scattering. 
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Fig. 6. — Source frame bolomctric luminosities for all 9 PI SNe out to 3 yr. Left panel: u-series. Right panel; z-series. The general trend 
of higher luminosity with progenitor mass in each series is evident, as is the fact that u-series explosions are brighter than z-series SNe of 
equal progenitor mass. The resurgence in luminosity at ~ 10^ s in most of the SNe coincides with the descent of the photosphere of the 
shock into the hot ^®Ni layer of the ejecta. 



several reasons for this: first, we envelop the star in a 
wind profile which is raised to extremely high temper- 
atures when the shock crashes through the surface of 
the star. This collision may account for some of the 
higher breakout l uminosity, since the re is no circumstel- 
lar shroud in the iKasen at al.l (|2011[ ) simulations. Sec- 
ond, we use 2T radiation transport instead of IT trans- 
port, and when radiation and matter are no longer forced 
to be in thermal equilibrium more photons may be emit- 
ted by the flow. We also use LANL OPLIB opaci- 
ties instead of the Lawrence Livermore OPA L opaci- 
ties (jlglesias fc Rogerslll996l : IRogers et ai]|1996l PI in the 

http;//rdc. llnl.gov 



Kasen et al. models, and this may too account for 
some of the disparity in peak luminosities. Shock break- 
out has been the subject of num erous analytical stud- 
ies (IColgatc.igfllMatzner fc McKee.1999 : Nakar fc Saril 
120101: IPiro et all 120101: iKatz et alTl2012D and numerical 



studies (lEnsman fc Burrowslll992l:TBlinnikov et alll200( ; 



iTominaga et all 120091: iTolstovl l2010t IKasen et al.H2011 ) 

in the past thirty years. One aspect of shock breakout 
that is unique to our study is that it is redshifted from 
1 - 2 hr in duration at z ~ 20 to a day or more in the 
observer frame. 

4.2. Intermediate / Late Stages of the SN 
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Fig. 7.— Spectral evolution of the z250 PI SN. Left panel: fireball spectra at 1.47 X 10'' s (black), 2.77 X 10"* s (blue), 1.15 X 10^ s 
(red), 6.12 X 10^ s (green), and 2.30 X 10^ s (purple). Right panel: spectra at 6.12 X 10*^ s (blue), 1.10 X lO'' s (red) and 2.08 X 10^ s 
(green), the onset, peak and decline of the bump in bolometric luminosity at ~ 4 months. Note the multiple P Cygni profiles of prominent 
lines (dips in flux blueward of the line and peaks in flux redward of the line), such as those at ~ 2.5 X 10** A, 4.0 X 10* A, 6.5 X lO** A, 
and 10^ A at 1.1 X lO'' s. 




We show bolometric light curves out to three years for 
the u-series and z-series SNe in the left and right panels of 
Figure [6j respectively. PI SN luminosities are powered 
at early times by the conversion of kinetic energy into 
thermal energy by the shock, so they are far brighter 
than Type la and II SNe at this stage because they have 
much higher explosion energies. At later times their lu- 
minosity comes mostly from radioactive decay, and they 
are much brighter than other SNe because they synthe- 
size more ^^Ni: up to 50 Mq compared to 0.4 - 0.8 Mq 
and < 0.3 Mq in Type la and Type II SNe, respectively. 
They are brighter for longer times (3 years instead of 3 - 



6 months for Type la and II SNe) because radiation dif- 
fusion timescales in their more massive ejecta are much 
longer: 

tdr^ K^Mej'^E--^. (8) 

Here, k is the average opacity of the ejecta, M^j is the 
mass of the ejecta, and E is the explosion energy. The 
luminosity generally rises with progenitor mass within 
each series because the explosion energy and ^®Ni mass 
increase with stellar mass in both red and blue stars. All 
five u-series explosions exhibit a slower, more protracted 
decay in luminosity out to 4 months than the z-series SNe 
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and are more than an order of magnitude brighter over 
this interval. These profiles are consistent with those of 
Type II-p SNe, whose progenitors are also thought to 
be red giants with extended envelopes. Their greater 
luminosities over this interval are mostly due to the fact 
that u-series explosions create more ^^Ni than z-serics 
SNe for progenitors of equal mass (see Table [T]) . 

Except for zl75, there is prominent bump in all the z- 
series luminosities at about 4 months that lasts for about 
a year. As we discuss in greater detail below, this re- 
brightening occurs when the photosphere of the shock 
uncovers the hot ^^Ni layer deep in the cjccta and its pho- 
tons escape into the IGM. The magnitude and duration 
of the bump is proportional to the ^^Ni mass, which ex- 
plains its absence in zl75 because it creates only a tenth 
of the 5^Ni formed in the z200 SN. This feature is also 
present in the u-series explosions but at first glance does 
not seem to be as prominent as in the z-series. Closer 
scrutiny reveals that it is several times brighter than in 
a z-series SN of equal mass, and that it is not as promi- 
nent because the luminosity before and after the bump 
are higher. The peaks are higher in the u-series SNe 
partly because they form more ^®Ni and partly because 
the photosphere extends deeper into the ^^Ni layer upon 
reaching it, where gas temperatures are higher. The radii 
of these layers are roughly the same when the photo- 
sphere reaches them in u-series and z-series explosions of 
equal mass. 

We show the evolution of the spectra for the z250 PI 
SN from shock breakout to 2 yr in the left panel of Fig- 
ure [T] Two physical processes govern the evolution of 
the spectrum over time. First, as the fireball expands 
it cools, and its spectral cutoff advances to longer wave- 
lengths over time. Second, the wind envelope that was 
ionized by the breakout pulse begins to recombine and 
absorb photons at the high energy end of the spectrum, 
as evidenced by the flux that is blanketed by lines at the 
short-wavelength limit of the spectrum. At later times 
flux at longer wavelengths slowly rises due to the expan- 
sion of the surface area of the photosphere. Over the 
wavelength scale of this plot the many thousands of lines 
captured by the LANL OPLIB opacities in our SPEC- 
TRUM calculation are blended together in the prominent 
jagged spectral features throughout the spectrum. Given 
how line blanketing by both the ejecta and the wind shear 
off the spectrum at short wavelengths, it is clear that the 
common practice of fitting blackbodies to light curves to 
approximate spectra overestimates flux at high energies 
from which many photons are eventually redshifted into 
the NIR in the observer frame. 

In the right panel of Figure [7] we show spectra at the 
beginning, peak and decline of the bump in bolometric 
luminosity at - lO'^ s for the z250 PI SN. At 1.10 x 10^ 
s, the emission peaks at ~ 1000 A and our SPECTRUM 
calculation reveals the r = 1 surface for this wavelength 
to be at the outer edge of the ^^Ni layer in the ejecta, as 
described above. We also note that these spectra mani- 
fest multiple P Cygni profiles for prominent lines, such as 
those at - 2.5 x 10"^ A, 4.0 x 10'' A, 6.5 x lO'' A, and 10^ 
A at 1.1 X 10^ s. They are characterized by the absorp- 
tion of flux blucward of the line and the enhancement 
of flux redward of the line. This occurs when line emis- 
sion from a fluid parcel receding along the line of sight 



from the observer is redshifted while line emission from a 
fluid parcel approaching the observer is blueshifted and 
absorbed en route to the observer. The presence of these 
features in our spectra validates our treatment of rela- 
tivistic time dflation and Doppler shifts in SPECTRUM. 
While we have singled out the z250 explosion for detailed 
spectral analysis, all the PI SNe spectra exhibit similar 
features and evolution. 

We plot velocity and density profiles at 5.02 xlO^ s, 
1.80 X 10'' s, and 1.0 x 10^ s for the u250 explosion in 
the left and right panels of Figure [51 Both forward and 
reverse shocks are visible in the velocity profiles until 
nearly the end of the simulation, after which the profiles 
become essentially self-similar. The reverse shock forms 
because the forward shock must propagate through an 
extended red hypergiant envelope before exiting the star, 
and it persists until well after shock breakout. Rayleigh- 
Taylor instabilities can arise in the reverse shock and lead 
to minor to modera t e mix ing in such stars, as sho wn by 
IJoggerst fc WhalenI (l20Tll) and IChen etall ([20T1 . Re- 
verse shocks also form in z-series SNe but they are not 
as strong because the stars are more compact and they 
have less time to develop before breakout. Their hydro- 
dynamical profiles become self-similar at earlier times. 

5. POP III PI SN DETECTION THRESHOLDS 

We have calculated NIR light curves for ou r PI SNe 
with t he synthetic photometry code described in lSu et al.l 
(|2011f ). Each spectrum is redshifted to the desired value 
before removing the flux absorbed by intervening neutral 
hydro gen along the li ne of sight according to the prescrip- 
tion of lMadaiir()1995D . We then dim the spectrum by the 
required cosmological factors. A variety of instrument 
fllters can be easily accommodated by our code, which 
linearly interpolates the least sampled data between the 
input spectrum and filter curve to match the other. It 
has additional capabilities such as reddening by dust that 
are not used here. 

5.1. NIR Light Curves 

At z > 7, Lyman absorption by the neutral IGM 
will absorb most flux at wavelengths longer than 1216 
A. However, in principle the flreball could be brighter 
slightly blueward of the Lyman limit rather than red- 
ward, even with I GM absorption (see Figure 3 in 
IWhalen et aHl2012bD . At every rcdshift for each explo- 
sion, we calculate the NIR signal in JWST NIRCam fil- 
ters above and below the Lyman limit to find the filter 
in which the SN is brightest. We find that for all the 
redshifts and PI SNe in our study the explosion is most 
luminous just redward of 1216 A in the source frame. We 
show light curves for our PI SNe at z = 10, 15, 20 and 
30 in Figures |9] and [101 At each redshift, light curves 
are plotted only for the filter in which the SN will be 
brightest. The detection limit in all four filters is AB 
magnitude 32. All five u-series explosions are visible for 
over 1000 days, even at z = 30. The z-series SNe are dim- 
mer but are all still visible out to z = 30 for 100 - 300 
days. They have lower bolometric and NIR luminosities 
(and dim sooner) because they create less '^^Ni than do 
the u-series. None of the light curves reach the detection 
threshold for 50 days, or 1.5 days in the source frame 
at z — 30, so shock breakout cannot be seen from Earth. 
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Fig. 9. — JWST NIRCam light curves for the 5 u-series PI SNe. Redshifts (z = 10, 15, 20, 30) are noted m the upper right corner of 
each panel. The wavelength of the optimum JWST filter at each redshift is noted on the y-axis labels and the times on the x-axes are in 
the observer frame. The dashed horizontal lines at AB mag 32 and 27 arc JWST and WFIRST detection limits, respectively. 



The X-rays and hard UV in the breakout transient are 
absorbed by the neutral IGM. 

As noted carher, as the fireball expands and cools its 
spectral peak steadily advances to longer wavelengths. 
This is manifest in the signal in the five NIR channels 
in Figure [TTJ The signal peaks at later times at longer 
wavelengths as the temperature of the fireball falls with 
time. The luminosity persists for longer times at larger 
wavelengths because the shock emits at these energies 
for a larger fraction of its cooling time than at higher 
energies. The expansion and cooling of the explosion, 
together with cosmological redshifting, also accounts for 
the shift of the NIR peak to later times and longer filter 
wavelengths with redshift in Figures [S] and 1101 For ex- 
ample, at z = 10 the u-scries light curves peak from 100 



- 300 days, and at z = 30 they peak from 200 - 600 days. 
As expected, the peak magnitude of each explosion rises 
with redshift. The JWST NIR photometry limit is mag- 
nitude 31 - 32 and magnitude 29 - 30 for spectrometry. 
All five u-series PI SNe will be bright enough for spec- 
troscopy out to z = 30 but it will only be possible out 
to z 10 for all four z-series explosions. This is impor- 
tant because resolving the order in which lines appear in 
spectra over time could provide a powerful probe of how 
heavy elements are mixed in the ejecta and thus of the 

explosion engine itself. 

As noted in iWhalen et all ()2012b[ ). the NIR flux 
evolves on timescales of ~ 1000 days in the u-series and 
exhibits much more variability than the bolometric flux 
in the observer frame (which evolves over 40 - 90 yr) be- 
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Fig. 10. — NIRCam light curves for the 4 z-series PI SNe. 

cause of the expansion and cooling of the firebah. Such 
variabihty is the key to discriminating these events from 
primitive galaxies, with which they otherwise overlap 
in color-color space. Their NIR flux rises much more 
quickly than it falls, so it is easiest to detect them in 
their earliest stages but they exhibit enough variation 
over survey times of 1 - 5 yr to be identified at later 
stages as well. If Pop III PI SNe arc found in the NIR, 
they will be thousands of times brighter than the halo 
or primitive galaxy that hosts them. Indeed, if a z ^ 15 
object exhibited any variation in luminosity in a survey, 
in all likelihood it would be a SN in a primeval galaxy. 

5.2. WFIRST, Euclid & WISH 

Pop III PI SNe could be found in large numbers in 
future all-sky NIR surveys such as Euclid, the Wide- 
field Imaging Surveyor for High-Redshift ( WISH) and 



WFIRST, whose target sensitivities at 2 fim are AB mag- 
nitudes 24, 27, and 27, respectively. In Figure[9]it is clear 
that WFIRST and WISH will detect u-series explosions 
out to z ~ 15 - 20. Our calculations indicate that even at 
z = 7, all nine PI SNe will be above magnitude 25, so Eu- 
clid will only detect them below this redshift (although 
higher z might be possible with spectrum stacking). It 
may be that the optimal redshift range for locating Pop 
III PI SNe is z ~ 15 - 20 because of Lyman- Werner 
(LW) UV radiative feedback. LW backgrounds from the 
first stars are thought to destroy H2 molecules and sup- 
press cooling in primordial halos, causing them to grow 
more massive before their interiors can self-shield from 
LW photons, form H2 and host primordial star forma- 
tion. Detailed numerical simulations show that the larger 
virial temperatures at the centers of such halos elevate 
cooling rates per H2 molecule there by 2 orders of mag- 
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nitude, leading to higher central collapse rate s that favor 
the formation of very ma ssive stars (O'She a fc NormanI 
l2007tlWise fc AbeH[2008B . LW backgrounds sufficient to 
postpone Pop III star formation likely did not arise until 
z ~ 20, delaying baryon cooling and collapse in halos un- 
til z '--^ 15 - 20. Global feedback may therefore enhance 
Pop III PI SNe rates at slightly lower redshifts that are 
well within the range of WFIRST and WISH. 

5.3. LSST/Pan-STARRS 

Could high-z PI SNe be detected in all-sky optical 
and NIR surveys by the Large Synoptic Survey Telescope 
(LSST) or the Panoramic Survey Telescope & Rapid 
Response System {Pan-STARRS)7 Unfortunately, above 
2 ~ 10 any source frame wavelength that could be red- 
shifted into the optical is extinguished by Lyman absorp- 
tion, but LSST has a target sensitivity of AB magnitude 
22 in the Y band (0.95 - 1.070 ^m) (which could be ex- 
tended to 25 with spectrum stacking) and Pan-STARRS 
has a detection limit of AB magnitude 26 in the Y band. 
As we show in Figure [T^ all five u-series explosions will 
be visible to Pan-STARRS and marginally detectable by 
LSST at z ~ 7. As we discuss in greater detail below, 
it has been speculated that very massive Pop III stars 
could form down to z ~ 6 in pockets of metal-free gas; if 
so, they could be found by these two observatories. 

6. POP 111 PI SN DETECTION RATES 

Although JWST will clearly be sensitive enough to de- 
tect z > 30 PI SNe, will it encounter such explosions over 
reasonable survey times given its narrow field of view? 
Their detection in a given survey critically depends on 
their event rates, which in turn are governed by primor- 
dial star formation rates (SFRs) and the Pop III IMF. 
Many physical processes regulate the Pop III SFR over 
cosmic time. Metals and UV feedback from early genera- 
tions of stars are especially important, since Pop III stars 
can only form in pristine gas and LW photons can destroy 
the II2 required for baryons in halos to cool and collapse 



into stars (e.g.. iHaiman et al.l 119971: I Glover fc Brandl 
[200lHMachacek et al.ll200lHO'Shea fc Normanll2008n . 

6.1. Semi- Analytical Estimates 

The original estimates of Pop III PI SN rates were 
based on simple halo mass distributions and cosmolog- 
ical parameters that lead to first star formation at ear- 
lier e pochs than do more improved parameters today 
(e.g.. iKomatsu et al.ll201ll ). They predict event rates 
that range from 0.1 - 1.5 sq. deg~^ yr~^ at z 25 
(|Wise fc Abelll200l to 0.2 a nd 4 sq. deg"^ yr^^ a t 
z ~ 25 and 15, respectively (jWeinmann fc Lilivl I2005D . 
These rates exclude clustering and radiative and mechan- 
ical feedback between halos (e.g., iMesinger et al.l [20061 : 
IMacIntvre eral|[200l I Whalen et al.ll2008al . 120101) but do 
consider global LW UV backgrounds that can delay star 
formation to slightly later epochs. 

6.2. Numerical Simulations 

Cosmological simulations that incorporate chemical 
and radiative feedback in varying degrees of detail have 
produced more realistic Pop III SFRs for the first bil 
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They basically agree on event rates at high red- 



shifts but differ somewhat from the original analytical 
estimates, mostly because they use more recent cosmo- 
logical parameters. By modeling the rise of LW back- 
grounds an d taking a simple ap proach to early metal 
enrichment. iHummel et al.l (j2012D find a cumulative rate 
of ~ 0.5 - 5 sq. deg-^ yr '^ for Pop III PI SNe at z 
> 5. IJohnson et all (pOll find PI SN rates of ~ 0.3 
sq. deg~^ yr~^ over the same redshift range, which are 
slightly lower because their models employ more sophis- 
ticated prescriptions for chemical and mechanical feed- 
back by SNe and assume a less top-heavy Pop III IMF. 
They adopted a PI SN progenitor ma ss range of 140 - 
260 Mq (e.g., iHeger fc WooslevI I2002D . but new mod- 
els have since extended this ra nge down to ~ 65 Mq 
(IChatzopoulos fc Wheeled I20TI . This new lower limit 
revises the PI SN rate in IJohnson et al.l (j2012l) upward 
by a factor of a few because of their choice of a Salpeter- 
like slope for the IMF. 

6.3. z= 15 - 20 and > 25 

The cumulative PI SN rates reported bv lHummel et al] 
()2012D imply lO'^ - lO'^ events per yr per JWST NIR- 
Cam field of view (10 arcmin^). Thus, JWST should be 
able to find 5 - 10 PI SNe over the lifetime of the mission, 
but greater numbers may be found if more time is ded- 
icated to the surveys. We note that Lyman absorption 
above z ~ 6, which we include in our stu dy, may slightly 
reduce the detection rates predicted by IHummel et al.l 
(|20T1 . All-sky surveys will find far greater numbers 
of PI SNe, albeit at somew hat lower redsh i fts. G iven 
the PI SN rates rep orted by IHummel et~all (|2012D and 
IJohnson et all (l20T2h . up to ~ 10^ PI SNe per year could 
be found at 15 < z < 20 by WFIRST and WISH. 

6.4. z< 10 

As noted earlier, PI SNe may also be found in all-sky 
optical and NIR surveys by LSST or Pan-STARRS but 
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time (days) 

Fig. 12. — Y-band (1.0 fim) light curves for all 5 u-serics PI 
SNc at z = 7. The z-series explosions are too dim to be detected 
by LSST or Pan-STARRS at this redshift. The dashed horizontal 
lines at AB mag 27 and 25 are Pan-STARRS and LSST Y-band de- 
tection limits, respectively (these thresholds assume spectral stack- 
ing). 

at lower redshifts [z < 10) because of ext inction by Ly- 
man absorption (sec also lFrost et alJ[2009l for a recent PI 
SN detection campaign in Spitzer data). Thus, if very 
massive Pop III stars form at 6 < z < 10, these sur- 
veys may discover their SNe if their rates are sufficiently 
high. Several recent numerical simulations suggest that 
such stars could form in isolated pockets of mctal-frce gas 
at z <6 (^o rnatore ct al. 2007; Trenti & Stia vell i 200l 
iTrenti et al.l 120091 ) . and such environ ments have now 
been d iscovered at even lower redshifts iFumagalli et al.l 
()20lH) . Including LW backgrounds and a simple pre- 
scription fo^mfitalenrkhment based on a numerical sim- 
ulation, iTrenti et al.l ()2009D estimate a cumulative Pop 
III PI SN rate of ~ lO'^ yr"! dcg-^ at 5 < z < 10, 
which implies an all-sky rate of up to ^ 10'^ yr~^. 

Other cosmological simulations that include chemi- 
cal enrichment and mechanical feedback by SNe have 
found lower Pod III SFRs fe.g.. iTornatore et al.l [20071 : 
iMaio et"ani2011l) that irn ply PI SN rates that are smaller 
than ITrenti et al.l ()2009[ ). However, newer models with 
better treatments of LW feedback and chemical enrich- 
ment now suggest that PI SN rates in this redshift 
range ar e higher by roughly an order of magnitude than 
those of 'Tornator e et alJ (I2007D and IMaio et al.l ()2011l ) 
dHummel ct al. 2 0121 : iJohnson et al.l I2012D. Siiice the 
Pop III SFR down to z 7 reported bv lWise et aII(l2012D 
is in g ood agreement with that found bv IJohnson et al.l 
(|2012D . their results are also consistent with a high PI 
SN rateE3 Overall, the results of numerical simulations 
indicate that the all-sky Pop III PI SN rate at 6 < z < 
10 may lie between ^ lO'' - lO"' yr~^, which bodes well 
for their detection by LSST and Pan-STARRS. 

^'^ See also lAhn et al.l l|2012l ). who find a similar evolution of the 
LW background due to massive Pop III stars at high-z. 



7. CONCLUSION 

We find that Pop III PI SNe will be visible in deep field 
surveys by JWST out to z > 30 and in all-sky surveys 
by WFIRST out to z - 15 - 20. They occur at rates 
that are sufficient to appear in deep field searches but 
will be discovered in much greater numbers in all-sky 
surveys that can be followed up by JWST and ground- 
based instruments. It may also be possible to find PI 
SNe at lower redshifts in current surveys, for example 
as Lyman break dropouts in the Hubble Space Telescope 
(HST) CANDELS survey or with the new Subaru Hyper 
Suprime Cam. Such strategies will be the focus of future 
studies. 

Our calculations emphasize the detection of Pop III PI 
SNe in the first few years of the explosion, but could their 
remnants be detected at later times by different means? 
iWhalen et all ()2008bD found that PI SNe in ionized halos 
eventually emit up to half of the original energy of the ex- 
plosion as H and He lines as the remnant sweeps up and 
shocks the relic H II region. Unfortunately, the luminos- 
ity of these lines is too low and redshifted to be directly 
detected. However, PI SNc can also deposit up to half 
of their energy into CMB photons by inverse Compton 
scattering, a nd could impose excess power on the CMB at 
small scales fOh ct al.l l200l iKitavama fc Yoshidal [20051: 
[Whalen et al. 2008lJ~ The resolution of current ground- 
based CMB telescopes such as the Atacama Cosmology 
Telescope and South Pole Telescope approaches that re- 
quired to directly image Sunyacv-Zcldovich (SZ) fluctua- 
tions from individual Pop HI PI SN remnants, so future 
observatories may detect them. 

The extreme NIR luminosities of primordial PI SNe 
could contribute to a NIR background excess, as 
has been suggested fo r Pop HI stars themselves (i.e. 
iKashlinskv et al.ll2005D . X-rays from gas plowed up by 
PI SN remnants, together with radiation from early black 
holes, would also build up an x-ray b ackground at high 
redshifts. Indeed, IJohnson fc Khochfa r (2011) have de- 
termined that x-rays from Pop HI SN remnants may 
have accounted for ~ 10% of the reionizing photon bud- 
get at early times. New calculations show that PI SNe 
will probably not appear at 21 cm because of their lower 
event rates at high redshift and because their remnants 
will not emit enough synchrotron radiation to be di- 
rcctl y detected by existing or future 21 cm observato- 
ries "(jMeiksin fc WhalenI [20T1 ) . The imprint of Pop HI 
PI SNe on the CMB and NIR backgrounds will be ad- 
dressed in future studies. 

If, as some numerical simulations and stellar archae- 
ology suggest, gas in primordial halos fragmented into 
multiple Pop HI stars that were tens of solar masses in- 
stead of hundreds, then CC SNe also occurred in the 
primeval universe. Such explosions would be similar 
in brightness to those in the local universe today be- 
cause their central engines mostly depend on the struc- 
ture of the inner 3-4 Mq of the star, which does not 
vary strongly with metallicitv (|Chieffi fc Limongil 120041 : 
IWooslev fc Hege^[2007t IWhalen fc Frvedl2012| ). Because 
they are 100 times dimmer than PI SNe and have softer 
spectra, it is not clear how much of their luminosity will 
survive Lyman absorption by the neutral IGM at high 
redshift. However, because a dozen or more such stars 
may form in the halo, CC SNe may be more plenti- 
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ful than PI SNe, which would enhance their likeUhood 
of detection. Furthermore, if the star ejects a massive 
sheU before exploding, the SN ejecta will Ught up the 
shcU in the U V upon crashing into it (super luminous 
Type Iln SNe - ISmith fc McCravlIM)?! : Ivan Marie et"an 
[20m ). Such events can have luminosities that can rival 
those of PI SNe, and they might be visible at higher 
redshifts tha n Type II explosions (iMoriva et al.l 120101 : 
iTanaka et aI1l2012t IMoriva et al.ll2012[) . We are now cal- 
culating light curves and sp ectra for 15 - 40 Mcr, Pop III 
Type II and Type Iln SNe (jWhalen et al.ll2012cl ial). 

Finally, strong gravitational lensing by massive inter- 
vening galaxies and clusters at 2; '--^ - 1 could boost flux 
from Pop III SNe, more than compensating for Lyman 
absorption and impr oving prospects for their detection 
(jRvdberg et al.|[2012f ). The probability that flux from a 
Pop III SN would be Icnscd in an all-sky survey and the 
degree of magnification both depend on the event rate 
at the given redshift, and may be fairly low. We have 
performed preliminary calculations that place the likeli- 
hood of lensing of 2 '-^ 20 objects at ~ 1 - 5% for flux 
boosts of 2 - 5. Much higher boosts (10 - 300) are possible 
near the edges of massive clusters but with much smaller 
search volumes and lower probabilities of encountering 
high-2 SNe. We continue to develop and refine Markov 
Chain Monte Carlo ray-tracing models of strong gravi- 
tational lensing of z ~ 20 events, the highest redshifts 
ever attempted, in order to assess its potential to reveal 
primeval SNe and galaxies. Although strong lensing is 
not necessary for detecting PI SNe it may be key to find- 
ing CC SNe in protogalaxies prior to reionization, less of 
whose flux survives Lyman absorption but whose higher 
event rates may favor the magnification of this flux. 

The detection of primordial SNe will directly probe 
the Pop III IMF for the first time and reveal the en- 



vironments in which they form. Their event rates will 
also trace the evolution of the first stellar populations. 
Pop 111 PI SN explosions at z ~ 10 - 15 will mark the 
positions of primeval g alaxies on the sky which might 
not otherw i se be found by J WST or TMT (jGreif et al.l 
[200l[2010l : IWise et al.ll2"0ll . Their discovery wiU open 
our first direct window on the era of first light. 
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